Reversible Band Gap Engineering in Carbon Nanotubes by Radial Deformation 
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We present a systematic analysis of the effect of radial deformation on the atomic and electronic 
structure of zigzag and armchair single wall carbon nanotubes using the first principle plane wave 
method. The nanotubes were deformed by applying a radial strain, which distorts the circular 
cross section to an elliptical one. The atomic structure of the nanotubes under this strain are fully 
optimized, and the electronic structure is calculated self-consistently to determine the response of 
individual bands to the radial deformation. The band gap of the insulating tube is closed and even- 
tually an insulator-metal transition sets in by the radial strain which is in the elastic range. Using 
this property a multiple quantum well structure with tunable and reversible electronic structure 
is formed on an individual nanotube and its band-lineup is determined from first-principles. The 
elastic energy due to the radial deformation and elastic constants are calculated and compared with 
classical theories. 

PACS numbers: 73.22.-f, 62.25.+g, 61.48.+C, 71.30.+h, 73.40.Lq 
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I. INTRODUCTION 

Modification of electronic properties of condensed sys- 
tems by an applied external pressure or strain in the elas- 
tic range have been subject of active study. However, in 
most of the cases, the changes one can induce by the 
elastic deformation is minute even negligible due to the 
rigidity of the crystals. On the other hand, the situa- 
tion is rather different for single wall carbc 
(SWNTs) owing to their tubular geometryM 
SWNTs are highly flexible and have a very large Young's 
modulus. They sustain remarkable elastic deformations, 
and it has been shown that the structure and elec- 
tronic propert^p-imfier.g^^p m by the these 
deformations. NyoMBlliyOBBEa Similarly, signifi- 
cant radial deformation of SWNTs can be realized in the 
elastic range, whereby the curvature is locally changed. 
This way, zones with higher and lower curvatures relative 
to the undeformed SWNT can be attained on the same 
circumference. Hence, one expects that radial deforma- 
tion can induce important modifications ill rtbg-el ectjcaru r 
and conduction properties of nanotubesEJ ootJnJcj'Ell 

Tight-binding calculations have indicated that a 
SWNT may undergo an insulateri-Hietal transition un- 
der a uniaxial or torsional strainMLi Multiprobe trans- 
port experimentsEj on individual SWNTs showed that 
the electronic structure can be modified by bending the 
tube, or by applying a circumferential deformation. Em- 
pirical Extended Huckel calculations!!^ predicted that the 
conductance of an armchair SWNT can be affected by the 
circumferential deformations and a band gap can develop 
on a metallic armchair SWNT upon twisting. The effect 
of the radial deformation and squei c^ing faav i e bogm , inves- 
tigated by using various methods.Eilt 
ever, in spite of these theoretical studies 



13, a system- 
atic analysis of the effect of the radial deformation on the 




electrical properties has not been carried out yet. 

The objective of this paper is to provide a better under- 
standing of the effect of radial deformation on the elec- 
tronic band structure and elastic properties of SWNTs, 
based on the extensive first-principle (ab-initio) total en- 
ergy and electronic structure calculations with fully op- 
timized structures. In the next section, a brief review of 
the the first principle pseudopotential plane wave method 
that we used will be given. The effect of the radial de- 
formation on the atomic structure will be discussed in 
Sec. Ill A. We discuss th e elas tic properties of SWNTs 
under radial strain in Sec. Ill B . We show that the calcu- 



lated elastic deformation energies as a function of radial 
strain can be described very w ell within the classical the- 
ory of elasticity. In Sec. Ill C, we discuss the effect of the 



radial deformation on the electronic properties. We find 
that zigzag nanotubes are metallized unde r radial defor- 
mation in the elastic range. In Sec. HID, this property 
is exploited to realize various quantum well structures on 
a single nanotube with tunable electronic properties. We 
applied two different radial deformations to two adjacent 
regions of a (8, 0) nanotube to generate band offsets at 
the interface, which in turn lead to multiple quantum 



well structures. Our conclusions are given in Sec. IV 
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II. METHODOLOGY 

The first principles total energy and electronic struc- 
ture calculations have been performed using the pseu- 
dopotential plane wave methodEa —within the general- 
ized gradient approximation (GGA)E3. Calculations have 
been carried out within periodically repeating supercell 
geometry because of the necessity of using the periodic 
boundary conditions with the plane wave method. We 
used a tetragonal supercell with lattice constants, a sc , 



b sc and c sc . The lattice constants, a sc and b sc , are cho- 
sen such that the interaction between nearest neighbor 
tubes is negligible (the minimum C-C distance between 
two nearest neighbor tubes is taken as 6.2 A). The lattice 
constant along the axis of the tube, c sc , is taken to be 
equal to the ID lattice parameter, c, of the tube. The 
tube axis is taken along the z-direction, and the circular 
cross-section lies in the (x, y)-plane. In the ID Bril- 
louin zone (BZ), the wave vector k z varies only along the 
z-axis. 

Plane waves up to an energy cutoff of 500 eV are used. 
With this energy cutoff aad using ultra soft pseudopo- 
tentials for carbon atomsc3, the total energy converges 
within 0.5 meV/atom. In addition to this, finite ba- 
sis set correctionsEll are also included. Owing to the 
very large lattice constants of the supercell, a sc and b sc , 
k-point sampling is done only along the-tube axis. The 
Monkhorst-Pack special k-point schemeE3 with with 0.02 
A -1 k-point spacing resulting 5 and 10 k-points within 
the irreducible BZ of the tetragonal supercell are used 
for (n, 0) and (n, n) tubes, respectively. 
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FIG. 1: Top view of the undeformed and deformed nan- 
otubes. The arrows indicate restoring forces on the fixed 
carbon atoms. The in-plane mirror symmetry can be bro- 
ken depending on the rotational orientation of the tube (see 
(d)). The corresponding point groups of nanotubes are also 
indicated. 



III. RADIALLY DEFORMED NANOTUBES 

A. Geometric Structure 

The radial deformation that is treated in this study 
is generated by applying uniaxial compressive stress a yy 
on a narrow strip on the surface of a SWNT. In practice 
such a deformation can be realized by pressing the tube 
between two rigid and flat surfaces. As a result, the 
radius is squeezed in the y-direction, while it is elongated 
along the x-direction, and hence the circular cross section 
is distorted to the elliptical one with major and minor 
axis a and b, respectively. A natural variable to describe 
the radial deformation is the magnitude of the applied 
strain along the two axes 
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where Rq is the radius of the undeformed (zero strain) 
nanotube. We note that the point group of the unde- 
formed nanotubes is D n h or D nc i for n even or odd, re- 
spectively. Under radial deformation described above, 
the point group becomes 6*2/1 or Dih- (see Fig. [|). How- 
ever, depending on the nanotube orientation around the 
tube axes the in-plane mirror symmetry can be broken. 
For the (6, 6) tube, we studied several different orienta- 
tions in order to investigate the effect of mirror symmetry 
on the band crossing at the Fermi level. Three different 
orientations with point groups Civ , C% and D% are shown 
in Fig. |l|d. 



For different values of strain e yy , we carried out full 
structural optimization under the constraint that the mi- 
nor axis was kept fixed at a preset value. The strains are 
in the elastic range, since the deformed tubes relax back 
to the undeformed state when the applied strain is re- 
moved. The structural relaxation is done in following 
steps: first, depending on rotational orientation of the 
SWNT, either a single bond or a carbon atom at both 
ends of the minor axis are pressed towards each other by 
(1 — e yy )Ro and are kept fixed. Then, under this con- 
straint, the coordinates of the remaining atoms and the 
lattice parameter of the tube c are optimized. At this 
step, some resultant forces remain on the fixed atom(s) 
with components opposite to the applied strain as well 
as perpendicular to it. In the second step, the fixed bond 
lengths are optimized together with all the internal co- 
ordinates of the atoms and c parameter. Eventually, in 
the final fully relaxed structure the only remaining force 
is the restoring force, opposite to the applied strain, on 
the fixed atoms. All other force components on these 
fixed carbon atoms and all the forces on the rest of car- 
bon atoms are optimized to be less than 0.01 eV/A. Fig- 
ure^ shows the cross-sectional view perpendicular to the 
tube axis of the fully optimized undistorted and distorted 
SWNTs as well as the restoring force vector. 

Figure |^a shows the pair distribution function in a de- 
formed and undeformed (7,0) SWNT. The first peak in 
Fig. [2|a corresponds to the first nearest neighbor distance, 
which is slightly broadened without a shift of the peak 
position with deformation. This indicates that the C-C 
bond lengths (« 1.41 A) are practically unaltered under 
the applied strain. Similarly, the second peak in Fig. |^a 
is also slightly broadened, indicating a small effect of the 
distortion on the second nearest neighbor distances. The 
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FIG. 2: (a) Pair distribution; (b) bond angle distribution 
functions for the (7, 0) SWNT. Solid line is for undeformed 
SWNT while the dotted line is for radially deformed one. 



effect of the radial deformation becomes apparent only 
for the third and further nearest neigbor distances. 

Given that first and second nearest neighbor distances 
did not change significantly with the radial deformation, 
the only remaining degrees of freedom is the bond an- 
gle as clearly seen from the angular distribution function 
shown in Fig. ||b. The main peak around 120° does not 
change with applied strain, but the other peak a few de- 
grees below the main peak for undeformed tube splits into 
j new peaks where j is the number of peaks in the radius 
distribution of zigzag (7,0), (8,0) and (9,0) SWNTs. On 
the other hand, for the armchair (6,6) SWNT, although 
the main peak is not changed with strain, the second peak 
is broadened by a few degrees. One direct consequences 
of this observation is that for zigzag tubes the lattice 
parameter c decreases very slightly with radial strain, 
whereas it is almost constant for the (6,6) SWNT. 

In summary, the radial deformation does not have a 
noticeable effect on the first and second nearest neighbor 
C-C distances but it induces significant changes in the 
bond angles. This observation is therefore important and 
has to be taken into account in tight binding studies of 
SWNTs with radial deformation. 
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FIG. 3: a) The strain component e xx — (Ro — a)/Ro along the 
major axis as a function of applied strain e yy — (Ro — b)/ Rq . 
The slope is the in-plane Poisson ratio, v\\ . (b) Variation 
of the elastic deformation energy per carbon atom, (c) The 
restoring force on fixed carbon atoms. For (8, 0) SWNT, the 
force is scaled by 0.5 since it is only on one carbon atom, while 
for the other tubes it is on two carbon atoms. 



B. Elasticity 

In order to describe the in-plane elasticity and defor- 
mation of the SWNTs, we use first-principle calculations 
of the elastic deformation energy, i.e. the amount of en- 
ergy stored in a SWNT as a result of radial deformation, 
and the classical theory of elasticity. The relation be- 
tween stress and strain is given by generalized Hooke's 
law, for the radial deformation described in the previous 
section 



Cxx — — Cne a 



C 



and 



F, 



.4 



(3) 



(4) 



where F y is the restoring force applied on the surface 
area A. C\\ and C\i are the in-plane elastic stiffness 
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constants. Assuming the validity of the Hooke's law, the 
strain energy becomes a quadratic function of strain as. 



Et(0) 



(5) 



The in-plane Poisson's ratio, uu, relates e xx and e yy , from 
Eq.| 
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The strain components are plotted in Fig. ^|a. As pre- 
sented in Table | z/ii decreases with increasing nanotube 
radius and is slightly smaller than 1.0. Eqs. ^ and [| can 
be cast in a simpler form by introducing v\\ : 
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where Ed is the elastic deformation energy obtained from 
the difference between the total energies of radially de- 
formed and undeformed SWNTs expressed in Eq. ^ (i.e 
Et{s xx , £yy) ~ Et{0)). At this point we examine how the 
stress and the elastic deformation energy calculated from 
first principles compare with the linear and quadratic 
forms in Eqs. [?] and || obtained from classical theory. To 
this end, we plot Ejj and the corresponding restoring 
forces, F y , as a function of e yy in Fig. |^b and c, respec- 
tively. Interestingly, the quadratic form obtained from 
classical theory fits very well to the elastic deformation 
energy calculated from the first-principles. Hooke's rela- 
tion, and hence clastic character of the deformations, per- 
sists up to e yy = 0.25. It is also noted that the SWNT be- 
comes stiffer as R decreases. The variation of the restor- 
ing forces is expected to be linear in the elastic range. 
The restoring forces in Fig. ||c are in overall agreement 
with this argument, except the deviations at certain data 
points due to uncertainties in the first-principle calcula- 
tions, which are amplified because the force is a deriva- 
tive quantity. Calculated elastic constants are listed in 
Table |. It is interesting to note that there are discrepan- 
cies in the theoretical results for Young's modulus, due 
to the assignment of the thickness, h, of the tube wallEl 
Two commonly used values are 3.4 A (based on graphite 
interlayer spacing) and 0.6 A (based on the n orbital ex- 
tent). The wall thickness h can be estimated from the 
present radial deformation data first by calculating the 
volume, O from Eq. ||. Then, h is solved by assuming 
that the tube is a slab with thickness h. From this anal- 
ysis, we found that h is radius dependent and it decreases 
from 0.88 A for (7,0) tube to 0.74 A for (6,6) tube. 



TABLE I: In-plane elastic constants of SWNT's. All elastic 
constants are in GPa except vu which is unitless. C e // = 
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FIG. 4: Energy band structures of undeformed (left) and ra- 
dially deformed (e yy = 0.23) (right) SWNTs along the V - Z- 
direction: a) and b) (7,0), c) and d) (8,0), e) and f) (9,0), g) 
and h) (6, 6). Solid line is the Fermi level. 



C. Electronic Structure 

We now discuss in detail the electronic structure of 
SWNTs under applied radial strain. The calculated band 
structures of undeformed and radially deformed zigzag 
(7,0), (8,0), (9,0) and armchair (6 6) SWNTs are pre- 
sented near the Fermi level in Fig. [|. The band gaps of 
zigzag tubes reduce with applied strain, and eventually 
vanish leading to metallization. Figure |5| summarizes the 
variation of band gap and density of states at the Fermi 
level, V(Ep), as a function of the applied strain. For 
(7, 0) and (8, 0) SWNT's the band gaps decrease mono- 




FIG. 5: The variation of the band gap, E g (a), and the density 
of states at the Fermi level T>(Ef) (b) as a function of applied 
strain e vv . 



tonically and the onset of an insulator-metal transition 
follows with the band closures occurring at different val- 
ues of strain. Upon metallization T>(Ep) increases with 
increasing strain. The behavior of the (9, 0) tube is, how- 
ever, different. Initially, the band gap increases with in- 
creasing strain, but then decreases with strain exceeding 
a certain threshold value and eventually diminishes. For 
all these zigzag SWNTs the band gap strongly depend 
on the magnitude of the deformation, and E g is closed at 
13 %, 22 % and 17 % strain for (7,0), (8,0) and (9,0) nan- 
otubes respectively. Whereas, the armchair (6, 6) SWNT, 
which is normally metallic, remains metallic with a slowly 
decreasing T>(Ep) cveafar significant radial deformation. 
Earlier Delaney et alEiLJ showed that the tt* -conduction 
and 7r- valence bands of a (10, 10) tube which normally 
cross at the Fermi level with quasi linear dispersion, open 
a pseudogap in the range of ~ 0.1 eV at certain direc- 
tions of the BZ perpendicular to the axis of the tube 
owing to tube-tube interactions in a rope. The open- 
ing of the gap is-caused by the broken mirror symmetry. 
Lammert et al£3 pointed out the gapping by squashing 
(20,20) and (36,0) metallic tubes, since circumferential 
regions are brought into close proximity. Uniaxial stress 
of a few kilobars can reversibly collapse a small radius 
tube inducing a 0.1 eV gap, while the collapsed lacge ra- 
dius tubes are stable. In the study of Park et alB-i, the 
bandgap of the (5, 5) tube were monotonically increasing 
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FIG. 6: The variation of energy eigenvalues of states near 
the band gap at the T-point of the BZ as a function of the 
applied strain. The shaded region is the valance band. The 
singlet state originating in the conduction band is indicated 
by squares. 



probably due to bilayer interactions, since the separation 
of the two nearest wall of the tube became comparable 
to the interlayer distance of graphite. 

In order to explain the band gap variation of (n, 0) 
tubes, the energies of a few bands near the band gap 
are plotted as a function of strain in Fig. [| The singlet 
7r*-state in the conduction band shifts downwards in en- 
ergy much faster than the other states do with increasing 
strain. This is due to the increasing curvature with in- 
creasing radial deformation. Since the singlet 7r* -state 
lies below the double degenerate 7r*-states for both (7, 0) 
and (8, 0) SWNT's, their band gaps are closed monotoni- 
cally with increasing e yy . On the other hand, for the (9, 0) 
SWNT this singlet 7r*-state is above the double degener- 
ate 7r*-states. The increase of the band gap at the initial 
stages of radial deformation is connected with relatively 
higher rate of downward shift of the double degenerate 
7r-valence band relative to the 7r*-conduction band un- 
der low strains. Once the singlet 7r*-band, which shows 
faster decrease with strain, crosses the doublet conduc- 
tion band and enters into the gap, the band gap begins 
to decrease with increasing strain. 

Finally, we examined the effect of the radial deforma- 
tion on the charge density. In Fig. ^ we show the charge 
density of states near the band edges. The effect of the 
deformation is remarkable on the singlet state; charge 
moves from the low curvature regions to the high curva- 
ture regions as the strain is increased. Significant charge 
rearrangements with radial deformation can modify the 
chemical activity of the surface of the SWNT relative 
to foreign atoms and molecules. Since a SWNT can 
sustain large elastic deformations, it allows significant 
charge rearrangements on its surface. Hence, this effect 
can be used to control-chemical reactivity of specific car- 
bon atoms in SWNTsEI. 
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FIG. 7: Charge density of the highest valence band state at 
T-point of a undeformed (a) and deformed (b) (9, 0) SWNT. 
The panels (c) and (d) show the the singlet "conduction" band 
states. 



D. Strain induced Quantum Structures 

It is clear from the above discussion that the band 
gap of an insulating SWNT can be modified, and even 
an insulator-metal transition can be induced by radial 
deformation in the clastic range. If the applied deforma- 
tion is not uniform but has different strength at different 
zones of the tube, it renders variable electronic structure 
along the tube axis. For example, each zone of an in- 
dividual SWNT undergoing different radial deformation 
attains a different band gap. Owing to the band off-sets 
at the junction, quantum structures can be engineered 
on an individual tube. 

Experimental and theoretical methods have been pro- 
posed in the past to determine the band offsets, and 
hence to reveal the band diagram perpetuating along the 
superlattice axistS However, there are some ambiguities 
in the case of SWNT superlattices (A n B m ), which are 
formed by periodically repeating undeformed A regions 
(of n unit cells) and radially deformed B regions (of m 
unit cells): (i) The alignment of the valence and con- 
duction bands of the A and B regions and the resulting 
band-lineup is a complex process involving charge trans- 
fer between A and B, and also modification of the crys- 
tal potential in the deformed region. In fact, the above 
analysis of radial deformation has shown that the va- 
lence and conduction band edges can be lowered when a 
zigzag SWNT is radially deformed. Therefore, a realistic 
treatment of the band alignment requires self-consistent 
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FIG. 8: The planar averaged crystal potential along the axis 
of the A(jBt> superlattice nanotube. Dotted vertical lines show 
the interfaces. Circular cross-section of the undeformed (8,0) 
nanotube in region A and elliptical cross-section of radially 
deformed nanotube (e yy — 0.16) in region B are shown as 
inset. The potential averaged over the original unit-cells of 
nanotube is shown by circles and dashed line. This average 
potential scaled by 25 for clarity. 



calculation of the crystal potential, (ii) Even if the band 
diagram were known, it is not obvious whether the Effec- 
tive Mass Approximation (EMA) is applicable for an in- 
dividual, nommiformly deformed SWNT. Therefore, in- 
stead of applying EMA to a ID real space band diagram 
or quantum well structure, one has to perform electronic 
structure calculations on the (A n B m ) supercell. How- 
ever, ab-initio calculations become tedious for large su- 
percell size owing to many involved carbon atoms. Ear- 
lier, we used a tight binding method and showed that 
states at the band edges are confined in either A or B 
regions of a superlattices (A 8 B S ), (A*J3i 2 ) and (A12B4) 
formed on an individual (7, 0) SWNTEI 

In this study, we extend our earlier calculations of 
quantum structures on SWNTs and present an ab-initio 
analysis of band lineup of the (AqBq) superlattice formed 
on an individual (8, 0) SWNT. Similar to the previous 
model, here the A region is left undeformed, but the B 
region is radially deformed by e yy = 0.16. This system 
consists of 384 carbon atoms in a supercell involving 12 
original unit cells of the (8,0) SWNT. We performed a 
partial structural optimization for the (A§B§) superlat- 
tice, since the full optimization is not tractable within 
a reasonable computation time. We first optimized sin- 
gle unit cells in the A and B regions corresponding to 
€ yy = and e yy = 0.16 respectively. Then, we connected 
A and B regions smoothly by one intermediate unit cell. 
Finally, fully self-consistent electronic band structure cal- 
culations were carried out on this structure. In Fig. |§| 
we show the planar (xy) averaged self-consistent poten- 



tial, V c (z) = J s V c (r)dxdy/S. Here S is the rcy-cross 
section of the supercell. The alignment of the valence 
band edges between undeformed A and deformed B are 
revealed by first integrating the planarly averaged poten- 
tial Vx = J) V c (z)dz/l at each region, (X — A or B), over 
a length of original unit cell (i.e. I = c) along the tube 
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axisE3. In Fig. g| Vx is shown for both regions of the 
supercell which constitutes the reference level for band- 
lineup. In the next step, we determine the energy of the 
valence band edge from the and Vg° calculated for 
two different, uniform (infinite) (8,0) SWNTs (one unde- 
formed, the other uniformly deformed with e yy — 0.16). 
These are Ey A and Ey B . It is assumed that in the nan- 
otube superlattice Ey A and Ey B are unaltered. The 
band-lineup of the valence band is calculated from the 
difference AE V = V A - V B + Ey° A - E^ B . For (8,0) 
SWNT, we find AE V ~ 180 meV; the valence band edge 
of B is lower than that of A indicating a staggered band 
lineup. This result clearly shows that by applying peri- 
odic radial deformation on an individual semiconducting 
SWNT one can generate a quantum structure, where the 
band gap in the direct space along the tube axis under- 
goes a periodic variation which is continuously tunable 
and reversible. 



IV. CONCLUSIONS 

In this work we present an extensive first principle 
analysis of the effect of radial deformation on the atomic 
structure, energetics and electronic structure of SWNTs. 
We find that the energy band structure and the variation 
of the gap with radius (or n) differs from what one derived 



from the zone folded band structure of graphene based 
on simple tight binding calculations. More interestingly, 
the response of the energy bands around the band gap 
to the applied radial deformation is different for differ- 
ent bands. Depending on the relative position of these 
bands, the band gap displays different behavior under the 
radial deformation. In general, the band gap is reduced 
and eventually closed to yield an insulator-metal tran- 
sition under the elastic radial deformation. The strong 
dependence of the band gap on the applied strain, its re- 
versible and continuously tunable behavior are exploited 
to form quantum well structure on an individual SWNT. 
A first-principle calculation of the alignment of the va- 
lence band is presented. The deformation energy and 
elastic constants under the radial deformation are calcu- 
lated. We find that the strain energy due to the radial 
deformation can be fitted very well to the quadratic ex- 
pressions obtained from the classical theory of elasticity 
within the Hooke's law. 
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